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Nonequilibrium properties of short-coherence-length s-wave superconductors are analyzed in the 
presence of extrinsic and intrinsic inhomogeneities. In general, the lowest-energy configurations 
of quasiparticle excitations are topological textures where quasiparticles segregate into antiphase 
domain walls between superconducting regions whose order-parameter phases differ by n. Antiphase 
domain walls can be probed by various experimental techniques, for example, by optical absorption 
and NMR. At zero temperature, quasiparticles seldom appear as self-trapped bag states. However, 
for low concentrations of quasiparticles, they may be stabilized in superconductors by extrinsic 
defects. 
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I. INTRODUCTION 

In general, studies of superconductors emphasize their 
equilibrium properties as probed by linear response. 
Equally important are conditions where the superconduc- 
tor is driven far from equilibrium. A nonequilibrium state 
may be achieved, for example, by photoexciting quasi- 
particles or injecting them into the superconduc- 
tor through a tunnel junction [§|5j|(|. These experiments 
have revealed a variety of interesting phenomena, rang- 
ing from first-order superconductor-metal transitions, to 
various instabilities and spatially inhomogeneous states 
with a laminar structure where either superconducting 
phases with distinct energy gaps or superconducting and 
normal phases coexist. 

At finite temperature, quasiparticle dynamics in the 
nonequilibrium state is dominated by scattering with 
phonons and decay with phonon emission. These pro- 
cesses are characterized by the scattering time r s and 
the lifetime r*. While these time scales are long enough 
compared to H/Aq so that the superconducting energy 
gap A is sharply defined 0, they are also typically of 
the same order of magnitude, r s ~ r* ||. As a con- 
sequence, it is not obvious that the quasiparticles will 
equilibrate to a meta-stable state even under steady-state 
conditions. In contrast, when the lifetime of quasiparti- 
cles is the longest time scale, quasiparticles will reach 
such a state making it possible for new phenomena to 
emerge. A meta-stable state is obtained, if the system 
has a symmetry group which makes it possible for the ex- 
cited and ground states to transform according to differ- 
ent irreducible representations so that the quasiparticle 
excitations cannot decay. It is clear that in the absence of 
spin-orbit coupling the only practically meaningful sym- 
metry group is spin-rotational symmetry, because in the 
superconducting phase broken gauge symmetry destroys 
the charge conservation. Quasiparticle excitations whose 



Ferromagnetic 
metal 




FIG. 1. Schematic description of a tunnel-junction exper- 
iment where spin-polarized quasiparticles are injected from 
a ferromagnetic metal through an insulating barrier to a su- 
perconductor. In the superconductor, injected quasiparticles 
have formed a laminar structure of intervening superconduct- 
ing domains separated by antiphase domain walls between 
regions where the neighboring order-parameter phases are 
shifted by tt. 

spins are aligned along the same direction may be ob- 
tained by using a ferromagnetic metal as a source ||. 

In this Note, we examine various meta-stable configu- 
rations of quasiparticles and their signatures that might 
develop when spin-polarized quasiparticles are excited 
in s-wave superconductors. Our most important find- 
ing is that superconductors are unstable against a for- 
mation of antiphase domain walls into which the quasi- 
particles localize and that the local structure and nonuni- 
form spin density makes these topological textures acces- 
sible to various experimental probes. In particular, they 
produce a distinctive optical absorption spectrum that 
may serve as a unique signature of their presence. In 
addition, any probe that is sensitive to a local magneti- 
zation would lend further support. For example, NMR 
and fiSK might be suitable for this purpose. Figure 1 il- 
lustrates schematically a tunnel-junction experiment for 
generating and detecting antiphase domain-wall textures. 
Similar experimental construction has been suggested to 
demonstrate that spin and charge are transported by 
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separate quasiparticle excitations in a superconductor 
. We also consider quasiparticle bags which are non- 
topological states of quasiparticles associated with a local 
suppression of the order parameter that may appear in 
the presence of defects when the quasiparticle density is 
small enough. 

Our work is partially motivated by the fact that only 
few studies exist on self-trapped quasiparticle states in 
superconductors ]llL|l2] and that either quasiparticle-bag 
or antiphase domain-wall excitations are usually consid- 
ered as a curiosity and often disregarded as unphysical 
|0,[l3). Our purpose is to address the question of their 
existence in mean-field approximation and to examine 
possible experimental implications by focusing on quasi- 
one and two-dimensional superconductors which are re- 
alized in wires and films. Analogous questions have been 
studied in the context of antiferromagnets, although no 
detailed predictions regarding superconductors have been 
made HQ. 

II. FORMALISM 

Our starting point in describing quasiparticle excita- 
tions in an s-wave superconductor is the lattice formula- 
tion of electrons hopping between nearest-neighbor sites 
and interacting via an effective two-particle interaction, 

H=-\W £ 4(R + r)i(R)-^ nj (R) 

(Rr)<r R<r 

+ [/]>> T (RH(R). (1) 

R 

Here, is the electron operator with spin a, (Rr) de- 

notes nearest-neighbor sites separated by r, W is the half 
bandwidth (on a square lattice), \x is the chemical poten- 
tial, and the operator 7V (r) = ?/>J.(r)^ CT (r) is the conduc- 
tion electron number density for spin a. The strength 
of the pairing interaction U (< 0) is assumed to be in- 
termediate so that the mean-field approximation gives a 
qualitatively reliable description of the superconducting 
ground state and the low-energy excitations. Specifically, 
consider a two-dimensional lattice model where electrons 
can interact with randomly distributed defects. These 
defects can either have a magnetic moment or be non- 
magnetic. The model is defined by the effective Hamil- 
tonian H c s = Hq + H lmp , where Hq describes a BCS 
superconductor |]l5| and -ff; mp is the contribution due to 
impurities. In the mean-field approximation, 

Ha = -jWJ2 * t (R + r)f 3 *(R)-/i^* t (R)f 3 *(R) 

(Rr) R 

-5>(R)*t (R)fl *( R)) ( 2 ) 

Rr 

where A(R) is the superconducting gap function and as- 
suming that the pairing of electrons occurs in the spin- 



singlet channel. The operator \&(r) = [^(r) V>j(r)] T is a 
two-component Gor'kov-Nambu spinor, f a (a = 1,2,3) 
are the Pauli matrices for particle-hole degrees of free- 
dom, and fo is the unit matrix. In a translationally in- 
variant system, the BCS Hamiltonian reduces to 

^o = ^* k (ekT 3 -A f 1 )* k , (3) 

k 

where Ao = A(R) and = (V'kt ^Ikj,) 7 - The fermion 
operators in real and momentum spaces are related by 
the unitary transformation, tp a (r) = N~ x / 2 ^ k i/we 1 * 1 ' 1- , 
where N is the number of sites in the system. For a 
square lattice with the nearest-neighbor hopping, the 
single-particle energy relative to the chemical potential 
in the normal state is ek = — ^VF(cos k x a + cos k y a) — fi; 
a is the lattice spacing. In a uniform s-wave superconduc- 
tor, the energy spectrum of bare quasiparticle excitations 
is Ek = + Aq. Allowing the excited quasiparticles 
to relax, the energy spectrum and the order parameter 
must be modified, as will be discussed below. 

The interaction between the conduction electrons and 
the impurities in the superconductor is given by the 
Hamiltonian 

ff imp =£[V(r)n(r)+JS(r). S (r)], (4) 

r 

where n(r) = Ya n a (r) and s(r) = \ Y av ?Aj(r)f (Tl/ ^(r) 
are the conduction electron number density and spin den- 
sity operators. In the case of point like impurities lo- 
cated at sites r„, the potential (scalar) and magnetic 
scattering terms have the forms V(r) = Y n V n 5rr„ and 
S(r) = Y, n SrAr„- Typically the distribution of impu- 
rities is assumed to be random whereas the magnitude 
of scalar and magnetic scattering are constant, V n = V 
and w = JS/2, where S — |S n |. For later emphasis, it 
is useful to introduce here a particle-hole transformation 
generated by the operator fi : 

(r) -> *'(r) - (-l) r f!*(r). (5) 

At half filling (/i = 0), the BCS Hamiltonian Ho on 
a square lattice is invariant under this transformation. 
Moreover, if the impurity moments are aligned along the 
same direction and there is no potential scattering, the 
impurity Hamiltonian H lmp will also be invariant un- 
der the same transformation. Potential scattering and 
randomly oriented impurity moments break particle-hole 
symmetry of this nature. 

Given that the pairing of electrons occurs in the spin- 
singlet state, the superconducting order parameter (am- 
plitude) can be expressed in the form 

F(R, r) = ±^r 2 U(Vv(R + r)Vv(R)). (6) 

(TV 

The relation between the order parameter and the gap 
function is given by the equation 
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A(R) = —UF(H,r = 0). (7) 

The on-site pairing interaction U is assumed to be in- 
stantaneous in time. Thus, the energy cutoff in the gap 
equation is set by the bandwidth. In our numerical ap- 
proach, the gap equation is solved self-consistently with 
a given number of quasiparticle excitations on finite size 
lattices with periodic boundary conditions. In our nu- 
merical examples, the strength of the interaction U is 
chosen so that in the absence of impurities and quasipar- 
ticle excitations the energy gap is A /W = 0.1. 

III. MAPPING TO AN ANTIFERROMAGNET 

The Hamiltonian (1) can be transformed to a model 
where the on-site interaction is repulsive. In the case of 
longer-range interactions, Ising-like terms are generated. 
On a bipartite lattice, this is achieved by a particle-hole 
transformation on the down spins, 

V> T (r) -> VtO), 
^(r) -> (-l)'^(r). 

In this transformation, the particle number operator 
transforms to the ^-component of the spin density opera- 
tor: n(r) — * 2s z (r) + 1, and vice versa. The Hamiltonian 
is mapped into 

H = -\W J2 V4(R + r)?MR) + ^$> 2 (r) 

(Rr)u r 

+ 2L/^,s z (r) Sz (r), (8) 

r 

where h z = U — 2fi is an effective magnetic field along the 
z-axis. Thus, for the Hubbard model, the particle-hole 
transformation changes the sign of the on-site interaction 
U. 

The superconductor has U(l) symmetry associated 
with the phase of the order parameter. Because the real 
and imaginary parts of the order parameter are trans- 
formed to the x and y components of the spin, a gauge 
transformation corresponds a rotation of the spin in the 
xy plane. 

The particle-hole transformation establishes one-to- 
one correspondence between the ground states of the at- 
tractive and repulsive Hubbard models. For example, 
consider the attractive Hubbard model away from half 
filling so that the average electron density (n) < 1 and the 
average spin density (s z ) = 0. The particle- hole trans- 
formation maps it into the half-filled, repulsive Hubbard 
model with the effective magnetic field h z . Its ground 
state has a transverse antiferromagnetic order because in 
this way the system can lower its energy by generating 
a small ferromagnetic component parallel to the z-axis. 



Therefore, in the ground state, (n) = 1 and (s z ) < 0. Re- 
versing the transformation, the transverse antiferromag- 
netic order parameter is mapped to a superconducting 
order parameter in the attractive Hubbard model. 

Next, consider the attractive Hubbard model away 
from half filling but now in the magnetic field so that 
(n) < 1 and (s 2 ) < 0. This model is mapped by the 
particle-hole transformation to the repulsive Hubbard 
model away from half filling. It has a ground state which 
is described by antiphase domain walls between antiferro- 
magnetically ordered spins [ p"3"[ . By virtue of the particle- 
hole transformation, it is clear then that the attractive 
Hubbard model with spin-polarized quasiparticles has a 
superconducting ground state where the superconducting 
domains with the opposite signs of the order parameter 
are separated by antiphase domain walls into which the 
excess spin is localized. 

IV. THE CONTINUUM MODEL 

Although we are mostly interested in quasi-two- 
dimensional superconductors, it is useful to consider one- 
dimensional systems where many ideas can be examined 
analytically. Indeed, for quasi-one-dimensional systems, 
a fruitful connection between the BCS and SSH Hamilto- 
nians can be made. The latter one describes, for example, 
conducting polymers where the order parameter A (a;) 
represents the lattice distortion |16|| . In our case, such 
systems can be organic superconductors or wires whose 
thickness is smaller than the coherence length £o- In the 
weak-coupling limit, additional progress is achieved by 
considering a continuum field theory, which can be de- 
rived because the coherence length is much longer than 
the lattice spacing, £o ^ a. 

Since we are interested in low-energy and long- 
wavelength phenomena, the electronic degrees of freedom 
can be expressed by slowly varying fields ip± a {x) describ- 
ing the left (+) and right (— ) moving electrons, 

Mx)/V^ = iP+AxV kFX + i>-*(x)e- ikpx , (9) 

where kp is the Fermi wave vector. Defining the four- 
component spinor as ^(x) = [$|(a;) $J(a;)] T , where 
= [ip +a (x) iii s (x)] T , the BCS Hamiltonian (2) 
becomes 

H Q = Jdx^(x)[v F p-A(x)n]^{x). (10) 

The momentum operator is p — —thf$d x , where vp = 
(2ta/h) sin kpa is the Fermi velocity. Similarly, the gap 
equation can be rewritten in the form p"7| 

A{x) = -\aU{^{x)f 1 ^{x)). (11) 

These equations are formally equivalent to those of the 
TLM model llq], which is the continuum limit of the 
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SSH model. For example, at zero temperature, the su- 
perconducting energy gap is A = 2IT / e~ 1 / A , where the 
dimensionless interaction is A = Np\U\; Np is the den- 
sity of states at the Fermi energy in the normal state. 
Similarly, the coherence length is £n — ^vf/Aq. 

It is now straightforward to determine the nonequi- 
librium properties of the quasi-one-dimensional s-wave 
superconductor. In particular, it is obvious that inject- 
ing spin-polarized electrons into the system, they form 
solitons. They are topological excitations of the system, 
acting as domain walls between two ground states that 
differ by the sign of the order parameter A. The en- 
ergy of the soliton is Ed w = 2A Q /n and the order pa- 
rameter is A(x) = Aotanh[(x — Xo)/£o]j where xq is the 
location of the center of the soliton. At low densities 
of injected electrons, single quasiparticle bags may ap- 
pear. They are counterparts of polarons; thus, also their 
spatial form as well as their energy, E qp = \/2Ed w , is 
known exactly. While in inhomogeneous superconductors 
individual quasiparticles may diffuse until they become 
trapped into defects, they are not generic solutions, be- 
cause at finite concentration of quasiparticle excitations 
they "phase separate" forming domain walls. 

V. ANTIPHASE DOMAIN WALLS 

While at zero temperature quasiparticle bags are not 
generic excitations of the superconductor, they may ap- 
pear as long-lived states because of defects. This may 
happen if they are injected into the system at low rate 
so that they can migrate without scattering from other 
quasiparticle excitations long distances before they are 
trapped to defects. Note that, in addition to magnetic 
impurities, a local order-parameter suppression caused 
by nonmagnetic impurities leads to bound states in the 
superconducting energy gap, albeit their binding energies 
are necessarily small JlJj] . Figure 2 illustrates a situation 
which is obtained when the quasiparticle concentration 
is small and there are magnetic impurities in the system. 
For simplicity, the magnetic impurities are assumed to 
be ferromagnetically ordered producing a maximal trap- 
ping potential. The bound quasiparticle states yield two 
peaks in the density of states, 

Af(w) = -^VlmG ffff (r,r;o; + iO+) ) (12) 

in the energy gap. The oscillations in Af(oj) for \u>\ > Ao 
are due to the finite size effects. 

With an increasing quasiparticle concentration, indi- 
vidual quasiparticle excitations become unstable towards 
a spontaneous formation of antiphase domain walls. This 
tendency is depicted in Fig. 3, where the number of quasi- 
particles is not large enough to form a domain wall that 
would extend all the way through the system. Instead, 
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FIG. 2. (a) The energy gap, (b) the spin density, and (c) 
the density of states of a localized solution of spin-polarized 
quasiparticles injected into an s-wave superconductor with 
magnetic impurities. This configuration of well separated 
quasiparticle bags is obtained self-consistently on a square 
lattice with the lattice spacing a, Aq/W = 0.1, ttNfw = 0.3, 
and \x = 0. The concentration of quasiparticles and magnetic 
impurities equals 1%. 

a closed domain-wall loop is formed. Because the order 
parameter changes sign across the domain wall, there are 
midgap states. The finite length of the domain wall leads 
to the level repulsion yielding the density of states that 
has a minimum at zero energy. As the system is half filled 
and either there are no impurities or their moments are 
parallel to each other, the effective Hamiltonian H c s is 
invariant under the particle- hole transformation, Eq. (|^). 
Consequently, the density of states, depicted in Figs. 2 
and 3, is symmetric relative the zero energy. 

For a finite concentration of quasiparticles, it becomes 
energetically favorable to form domain walls with infinite 
length; see Fig. 4. This allows all the quasiparticles to 
occupy the midgap states. Domain walls may become 
pinned to defects either because the defects have a mag- 
netic moment or because the defects suppress the order 
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FIG. 3. (a) The energy gap, (b) the spin density, and 
(c) the density of states when a finite number (40) of 
spin-polarized quasiparticles is injected into an s-wave super- 
conductor. This configuration is obtained self-consistently on 
a square lattice with the lattice spacing a, Aq/W = 0.1, and 
(i — 0. No impurities are present. 

parameter locally, and this local suppression then pins 
a domain wall. In the case of extended defects, domain 
walls may find it preferable to wind through these de- 
fects. 

All the quasiparticle and domain-wall textures are 
charge neutral when the system has particle-hole sym- 
metry at the Fermi energy. In this regard, quasiparti- 
cle bags can be described as spinons [p~of] , because they 
carry spin but no charge. On a square lattice with the 
nearest-neighbor hopping, this happens exactly at half 
filling (/i = 0). However, if particle-hole symmetry at 
the Fermi energy is broken, self-consistently determined 
quasiparticle configurations usually acquire charge, be- 
cause they are a linear combination of plane-wave states 
with an energy spread Ae ~ ?ivf/£o about the Fermi 
energy. Similarly, away from half filling, the domain 
walls become charged, although their total charge per 
unit length can be quite small. This feature is naturally 
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FIG. 4. (a) The energy gap, (b) the spin density, and 
(c) the density of states in an s-wave superconductor with 
5% spin-polarized quasiparticles. The antiphase-domain-wall 
configuration is obtained self-consistently on a square lattice 
with the lattice spacing a, Ao/W = 0.1, /i = 0, irNpV = 0.3, 
and rii m p = 2%. 

understood by considering the repulsive Hubbard model 
with a finite effective magnetic field which induces a small 
longitudinal ferromagnetic component. In the supercon- 
ductor, this component is equivalent to a non-zero charge 
density. In contrast, bare quasiparticle excitations at the 
Fermi surface (k = k^) behave as spinons irrespective of 
the energy spectrum in the normal state. 

Finally, consider the stability of domain- wall solutions 
against a formation of isolated quasiparticle bags. Their 
energies per particle can be computed numerically. In 
two dimensions, the energy of a vertical domain wall 
per particle is estimated as E^ w ~ 0.66A and the en- 
ergy of a quasiparticle bag as E qp ~ 0.86A . These es- 
timates are in agreement with those computed in the 
antiferromagnetic system for vertical domain walls 
and spin polarons po|| . Thus, approximately at the tem- 
perature T* ~ An/5 a considerable fraction of domain 
walls begins to evaporate forming isolated quasiparticle 
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bags. It is interesting to compare this temperature with 
the critical temperature of the superconductor, which is 
T c ~ Ao/2. Thus, there is a sizable temperature regime 
below T c where most of the excitations appear as iso- 
lated quasiparticle bags. At low enough temperatures, 
T < T c /3, domain-wall textures are thermodynamically 
favored over non-topological quasiparticle configurations. 



VI. OPTICAL ABSORPTION 

The optical absorption provides a specific probe to var- 
ious inhomogeneous states of nonequilibrium supercon- 
ductors. The optical absorption is the real part of the 
complex conductivity, a' ab (u>) = Kea a b(u>) (a, b = x,y), 
where 

1 nc^ 
a ab {uo) = — [A ab (q = 0, uj + i0+) + — Sab] ■ (13) 
iui 1 m* 

The ratio between the density of charge carriers and their 



effective mass is defined as n/m* 



(H kin )/2, where 



ifkin is the kinetic-energy part of the Hamiltonian H. 
The current-current correlation function is given by the 
formula 



A o6 (q,i) = -(Tj a (q,f)j fe (-q,0)), 
and its Fourier transform is 



(14) 



1 f°° 

A ab (q,u) = -J dte l ^A ab {q,t). (15) 

The current operator in Heisenberg picture is defined as 
j a (q,£) = e lHt j a (q)e- lHt : where 

j a (q) = \eaW ]T *t( R + r a )f *(R) e ^ r . (16) 
<Rr) 

It is useful note that optical absorption obeys the sum 
rule: 



7r ne 
2 m* 



(17) 



It is a quantity describing any state that is linearly per- 
turbed by the electric field. Typically, it is associated 
with the equilibrium state. For non-equilibria states, 
such as the domain walls and quasiparticle bags, the sum 
rule must be modified. 

In the normal state, the optical conductivity has the 
Drude form 



ne t 1 

TO* (tw) 2 + 1 ' 



(18) 



where r _1 is the scattering rate due to the impurities. In 
the limit of dilute concentration of impurities, it can be 
approximated as 
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FIG. 5. Optical absorption a' xx (to) in an s-wave supercon- 
ductor in the presence of quasiparticle excitations forming 
a domain-wall lattice and localized quasiparticle bags. The 
concentration of quasiparticles in both cases is 1%. These 
configurations are obtained self-consistently in one dimension 
for A /W = 0.1, n = 0, nN F V = 0.3, and n imp = 5%. 
The dashed line denotes the optical absorption (Drude-like) 
obtained in the normal state (Aq = 0). 
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(19) 



where <5 is the phase shift for s-wave scattering and ni mp 
is the impurity concentration. For point-like impurities, 
the scattering phase shift is obtained from the equation 
cot 5 = c, where c = (irNpV^ 1 , for nonmagnetic im- 
purities, and c = (irNpw)^ 1 , for magnetic impurities. 
Below, as a reference, the numerically determined opti- 
cal conductivity in the presence of randomly distributed 
impurities in the normal state (Ao = 0) is also shown. It 
is well described ^ij by the Drude form, Eq. (p"8|). 

In the superconductor, the zero-temperature optical 
conductivity typically has a threshold of 2Ao due to the 
superconducting energy gap in the electronic spectrum at 
the Fermi energy. Quasiparticle bags and antiphase do- 
main walls introduce states within this energy gap that 
can be used as a characteristic signature of them. Fig- 
ure 5 illustrates the optical absorption when the injected 
quasiparticles form either isolated bag states pinned to 
nonmagnetic impurities or domain walls. In the former 
case, there is a very large peak in the absorption that 
comes from the excitation processes from the localized 
bag states to states just below the energy gap. Because 
the order-parameter suppression occurs on the length 
scale determined by the coherence length, it acts as an at- 
tractive potential with a finite range that can bind states 
just below the energy gap. Thus, in addition to the state 
occupied by the quasiparticle, the order-parameter re- 
laxation may admit additional discrete states below the 
energy gap. Transitions between these states have a very 
large oscillator strength. One can also see a peak at 2Ao, 
which is due to the pair-breaking processes across the 
energy gap. In the case of domain walls, the optical ab- 
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FIG. 6. Optical absorption o' xx (lo) in an s-wave super- 
conductor with 5% scalar (ttN f V = 0.6) and magnetic 
(ttNfw — 0.6) impurities in the absence of quasiparti- 
cle excitations. The ground-state configuration is obtained 
self-consistently in one dimension for Ao/W — 0.1 and /i = 0. 
The dashed line denotes the optical absorption (Drude-like) 
obtained in the normal state (An = 0). 

sorption begins at An due to the midgap states. One can 
therefore clearly distinguish between these two nonequi- 
librium states based on the optical absorption. 

Impurities yield a quite different absorption spectrum 
in the absence of quasiparticlc excitations; see Fig. 6. 
Nonmagnetic impurities produce a spectrum that has a 
clear threshold near 2A and a peak, whereas magnetic 
impurities yield a relatively smooth absorption profile 
which may extend deep below 2A , depending on the 
coupling strength between electrons and impurity mo- 
ments. For a high enough concentration of magnetic 
impurities, the superconductor becomes gapless and the 
absorption will begin at zero energy. 



This will require the use of spin-polarized quasiparticlcs, 
which may not always be feasible. A qualitatively similar 
situation may be created by maintaining a steady state 
of unpolarized quasiparticlcs by continuously pumping 
quasiparticles into excited states. Even though a gen- 
uinely metastable state may not develop because r* ~ r s , 
the fact that t* can be many orders of magnitude longer 
than the time scale associated with the superconducting 
energy gap, h/A , suggests that some of the features ex- 
plored here may actually be relevant for such states, too. 
Time resolved techniques are an ideal tool to probe their 
properties. 

While in gapless superconductors, such as in d-wave 
superconductors, it is no longer clear that quasiparticle 
excitations will lead to antiphase domain walls, various 
external defects that suppress the order parameter may 
locally favor a phase shift. Such textures may appear in 
magnetic superconductors with static spin-density-wave 
ordering where the phases of the magnetic and super- 
conducting order parameters intertwine to form a new 
collective state with midgap quasiparticle states. 
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VII. FINAL REMARKS 

Based on both analytical and numerical approaches, 
we have demonstrated that s-wave superconductors 
driven away from equilibrium exhibit interesting topo- 
logical textures. They develop as quasiparticles in node- 
less superconductors segregate forming antiphase domain 
walls in the superconducting order parameter and in this 
manner induce low-energy excitations into which quasi- 
particles relax. Their inhomogeneous structure has clear 
experimental implications. For example, a nonuniform 
spin density associated with domain walls should be ac- 
cessible to any probe that is sensitive to a spatially vary- 
ing magnetization. Moreover, optical absorption pro- 
vides another unambiguous tool for exploring these tex- 
tures. 

We have assumed that the lifetime r* of the quasi- 
particles in the excited state is much longer than the 
scattering time r s so that a metastable state is reached. 
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